Fitting inverse power-law quintessence models using the SNAP satellite 
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We investigate the possibility of using the proposed SNAP satellite in combination with low-z 
supernova searches to distinguish between different inverse power-law quintessence models. If the 
true model is that of a cosmological constant, we determine the prospects of ruling out the inverse 
power-law potential. We show that SNAP combined with e.g. the SNfactory and an independent 
measurement of the mass energy density to 17% accuracy can distinguish between an inverse power- 
law potential and a cosmological constant and put severe constraints on the power-law exponent. 
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I. INTRODUCTION 

Recent high precision measurements on type la super- 
novae (SNIa), the cosmic microwave background (CMB) 
and rich galaxy clusters indicate that our Universe is well 
described by a flat Friedmann-Robertson- Walker (FRW) 
model with mass energy density Q m ~ 0.3 and dark en- 
ergy density fid ~ 0.7 M. However, our knowledge of 
the nature of the dark energy is still limited, amount- 
ing to the facts that it is smooth on small scales and 
has negative pressure, therefore causing the Universe to 
accelerate. The most famous dark energy model which 
meets these two requirements is vacuum energy due to 
the cosmological constant A, but equally well-known are 
the many problems associated with its introduction 
A fairly recent alternative to A is quintessence, which re- 
sembles inflation in that it introduces a new scalar field 
Q to account for the dark energy. Consequently, the 
equation-of-state parameter wq = pq/pq is in general 
time- varying, whereas wa = — 1 always. At the present 
epoch, the field sits still or rolls slowly down its potential 
which results in a negative equation of state, i.e. nega- 
tive pressure. A key feature of many quintessence models 
is that the equations of motion have attractor-like solu- 
tions, which ensures that a wide range of initial condi- 
tions give the same late-time evolution of the field. Al- 
though this behaviour solves one of the fine-tuning prob- 
lems of the cosmological constant, the tuning of initial 
conditions, it does not answer the question why fi m ~ fid 
today @. 

The prospects of using the proposed Super- 
Nova/ Acceleration Probe (SNAP) [Q, a two-meter 
satellite telescope dedicated to the search and follow- 
up of supernovae, as a tool to probe the nature of 
dark energy has been examined in several papers (see 
@, | @, | § @ and references therein). For 
example, Goliath et al. show that given an inde- 
pendent high- precision measurement of Q m , SNAP 
can constrain the parameters of a linear equation of 



state wq{z) = Wq + WqZ to within ±0.04 and t^rr 
respectively at the one-sigma level, assuming a flat 
Universe. Moreover, Weller & Albrecht conclude that 
SNAP will be able to distinguish some quintessence 
models from a pure cosmological constant. In this paper, 
we investigate the possibility to differentiate between 
different parameters in one and the same potential. As 
our data we will use simulated SNIa magnitude-rcdshift 
measurements corresponding to one year of SNAP data 
and 300 events from low-z supernova searches such as 
the SNfactory @. 

One of the most common quintessence models in the 
literature is the simple inverse power-law potential in- 
troduced by Ratra and Peebles Q and reanalysed by 
Steinhardt et al. IU5|, 



V = 



M 4+c 



(1) 



Here M is a mass parameter which is fine-tuned to give 
the right fid today when Q is of the order of unity in 
Planck units 15 . The exponent a is the parameter we 
want to estimate. It determines the value of wq to- 
day, with smaller a giving larger negative wq. In the 
limit a — ► we retrieve the cosmological constant with 
PA = M A . Although there is some theoretical motivation 
for such a potential from supersymmetric QCD (see (ll| 
and references therein), the main reason for its introduc- 
tion is phenomenological. In addition to the attractor- 
like solutions mentioned above, potentials of this type 
produce an equation-of-state parameter wq which au- 
tomatically decreases to a negative value at the onset 
of matter domination. Since the energy density of any 
component evolves as pi oc a~ 3 ( 1 + , " i ) (a being the scale 
factor), this means that while quintessence starts out as 
a subdominant contribution, it will eventually come to 
dominate the Universe. 



II. METHOD 



A. Basic equations 
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We assume a spatially flat FRW Universe with zero 
cosmological constant and set the Planck mass to unity. 
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For a homogeneous and minimally coupled scalar field, 
the equations of motion are then given by 



Q + 3HQ + V'{Q) = 0, 
H 2 = ^-(PQ+ Pb), 



(2a) 
(2b) 



combined experiments. The individual measurement er- 
ror is assumed to be a — 0.1 5 mag, including the intrinsic 
spread of supernova magnitudes but ignoring systematic 
effects from e.g. gravitational lensing or dust. 



where H is the Hubble parameter, is the background 
energy density and dot and prime represent differenti- 
ation with respect to time and Q respectively. These 
equations can be solved numerically for the redshift de- 
pendence of the equation of state parameter 



V 



V 



(3) 



independently of the value of the Hubble constant Hq. 

The apparent magnitude m of a supernova is in- 
dependent through the luminosity distance d^, 



M + 5log(H d L ) 



(4) 



Here M. is a constant which depends on Hq and the ab- 
solute magnitude of the supernova which we assume not 
to evolve with redshift. The second term however, is Hq- 
independent and is given by 

H d L (z) = (1 +z) x 

x f dz' [n m {l + z'f +n Q f(z')]^ , (5) 
Jo 

,l + w(z') 



(6) 



At low redshifts, H^dL ~ z for any cosmology which 
means that low-z events can be used to measure M. with- 
out any prior knowledge of the cosmological parameters 
in general and Hq in particular. This is why Ai is some- 
times referred to as the Hubble-constant-free absolute 
magnitude. 



B. Simulating supernovae 

We use the SuperNova Observation Calculator 
(SNOC) @, a program package developed at Stock- 
holm University, to make Monte Carlo simulations of 
magnitude-redshift datasets attainable with SNAP and 
other experiments. The fiducial FRW cosmology we use 
is parametrized by M. = —3.4, Q m = 0.3, 17q = 1 — fl m = 
0.7 and the inverse power-law exponent a an integer num- 
ber in the range — 10 (with the value corresponding 
to a cosmological constant). 

SNAP is expected to observe ~ 2000 SNla per year 
for three years at redshifts out to z = 1.7 [Q. To this 
sample we add 300 events at z < 0.15 from future low-z 
supernova searches such as the SNfactory. This is proba- 
bly a conservative estimate of the amount of low redshift 
data available at the time of the SNAP launch in 2009. 
Figure [l] shows our simulated SNla distribution for the 
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FIG. 1: Simulated SNla distribution for one year of SNAP data 
combined with 300 low-z events from e.g. the SNfactory. 



C. Fitting procedure 



We want to fit our simulated magnitude-redshift mea- 
surements for a using the method of maximum likeli- 
hood. Since we have assumed a flat Universe, the other 
unknown parameters are f2 m and M. . The best- fit values 
are then found by minimizing the negative log-likelihood 
function C: 



1 n 

IE 



[m(6, Zi) - m(0 sim , Zi)Y 



(7) 



where 6 = (a,f2 m , A4) and subscript sim means simu- 
lated values. Note that additive terms which only result 
in a shift in C have been dropped in the above expression. 

The parameter ranges we consider are < a < 10.5, 
0.1 < fl m < 0.6 and —10 < M. < 0. In order to minimize 
C numerically, we need to know m and thus wq for any 
given 0. Therefore we first solve the equations of motion 
(||) for the evolution of wq using 165 different (a, f2 m ) 
combinations in the above ranges and then use linear 
interpolation to find wq for any combination of a and 
ft m . We will display our results using confidence regions 
in the (a, f2 m ) plane. 
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III. CONFIDENCE REGIONS 
A. No priors 

In this section we assume no prior knowledge of M. nor 
fi m . This means that for each grid point in the (a, Q m ) 
plane, we find the value of M. which minimizes the neg- 
ative log-likelihood function. Naturally, since we are not 
trying to estimate M, this value is of no importance to 
us. 

Figure || shows 68% confidence regions with one year of 
SNAP data. Different contours correspond to the fidu- 
cial cosmology with < a < 10. Due to the strong 
correlation between a and f2 m , it is practically impos- 
sible to distinguish between different integer a models. 
The plotted confidence regions are forced to be centered 
on the true values of the fiducial cosmology, but that will 
most certainly not be the case for the real confidence re- 
gion from only one experiment. However, the size and 
shape of the region will remain approximately the same 
which means that the region for any a ^ may drift 
to cover almost the entire plane. For a < 2, we can still 
estimate a to better than at the one-parameter, onc- 
sigma level. In the case of a cosmological constant and a 
measured value centered on a = 0, we can conclude that 
a < 0.9 at the one-parameter, one-sigma level. 
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FIG. 2: 68% confidence contours for (a, Q m ) with one year of 
SNAP data. The different contours correspond to the simulated 
fiducial cosmology with < a < 10. 

Combining SNAP with e.g. the SNfactory reduces the 
error contours somewhat as shown in figure |3j. The reason 
for this was given above; adding low-z supernovae corre- 
sponds to prior knowledge of M., something which will 
become evident later when we assume exact knowledge of 
M.. For a < 2, we can now estimate a to better than ^\ 
at the one-parameter, one-sigma level. In the case of a 
cosmological constant and a measured value centered on 
a = 0, we find a < 0.6 at the one-parameter, one-sigma 



level. However, to break the degeneracy in parameter 
space we need an independent assesment of O m . 
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FIG. 3: 68% confidence contours for (a, f2 m ) with one year of 
SNAP data and 300 low-2 events. The different contours corre- 
spond to the simulated fiducial cosmology with < a < 10. 



B. Priors on Q m and M 

Our assumption of no prior knowledge of Q m is actu- 
ally unrealistic, especially since we already assume that 
the Universe is flat, thereby leaning heavily on the latest 
results from CMB measurements. These same results, 
as well as independent observations of the evolution of 
the number density of galaxy rich clusters and mass esti- 
mates of galaxy clusters, all favour a low mass Universe. 
However, we must be certain that these measurements 
really are independent of wq before using them to in- 
voke prior knowledge of fl m . In plf , Gerke & Efstathiou 
give a brief summary of present and future experiments 
which may resolve this issue and conclude that at the 
time of the SNAP launch, wq -independent measurements 
of fi m corresponding to a Gaussian prior knowledge with 
a spread of <7n m = 0.05 may actually be a conservative 
expectation, whereas <Jn m — 0.015 could be considered 
optimistic. We note that a recent analysis of 2dF data 
already gives ft m = 0.27 ± 0.06 fT|. 

The result of imposing these priors respectively is 
shown in figures [| and ||. With crn m = 0.05, we can 
determine a to better than j^'g at the one-parameter, 
one-sigma level for all a. In the case of a cosmological 
constant and a measured value centered on a = 0, we 
can conclude that a < 1.2 at the one-parameter, three- 
sigma level. If we instead use the optimistic spread of 
fn m = 0.015, we can determine a to better than at 
the one-parameter, one-sigma level for all a and in the 
case of a cosmological constant and a measured value 
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centered on a = 0, we can conclude that a < 0.6 at the 
one-parameter, three-sigma level. 
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FIG. 4: 68% confidence contours for (a, Q m ) with one year of 
SNAP data and 300 low-z events, assuming a prior knowledge with 
H m Gaussian centered around 0.3 and <JQ m = 0.05. The different 
contours correspond to the simulated fiducial cosmology with < 
a < 10. 



be launched before 2012, is expected to find roughly 10 5 
low-z supernovae of all types within 4 years []l9f . 

Figure ^| shows 68%, 95% and 99% confidence re- 
gions for (a, f2 m ) assuming exact knowledge of M and 
&n m = 0.015. It is clear that such a scenario would en- 
able pinpoint precision in the determination of a and in 
the case of a cosmological constant the ability to rule out 
a > 1 at the 99% confidence level. 
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FIG. 5: 68% confidence contours for (a, Q m ) with one year of 
SNAP data and 300 low-z events, assuming a prior knowledge with 
Vim Gaussian centered around 0.3 and (Tq = 0.015. The different 
contours correspond to the simulated fiducial cosmology with < 
a < 10. 

Finally, we have performed the same analysis assum- 
ing exact knowledge of M. . This scenario is admittedly 
Utopian but it serves to illustrate the full potential of the 
SNAP satellite when combined with other experiments. 
Moreover, there are plans for future experiments which 
could provide values of M with a extraordinary preci- 
sion. The ESA satellite Gaia for example, planned to 



FIG. 6: 68%, 95% and 99% confidence contours for (a, Q m ) with 
one year of SNAP data assuming exact knowledge of Ai and a 
prior knowledge with Q m Gaussian centered around 0.3 and <JQ m = 
0.015. The different contours correspond to the simulated fiducial 
cosmology with < a < 10. 



IV. CONCLUSIONS 

We have investigated the possibility of using future 
SNIa observations to constrain the exponent a of the in- 
verse power-law quintessence potential. We do not use 
any parametrization of the equation of state nor any of 
the "standard practices" which according to |ll| may 
result in substantial estimation errors. In spite of the 
outstanding redshift range and measurement precision of 
the SNAP satellite, as well as the relatively accurate de- 
termination of M. possible with e.g. the SNfactory, it is 
clear that an independent measurement of the mass en- 
ergy density O m is needed in order to make any firm 
predictions of the value of a. 

However, armed with a Gaussian prior on f2 m with a 
spread of an m = 0.05, we can determine a to better than 
t\~3 at the one-parameter, one-sigma level for 1 < a < 
10. In the case of a cosmological constant and a measured 
value centered on a = 0, we can conclude that a < 1.2 
at the one-parameter, one-sigma level. In order to take 
full advantage of the SNAP satellite's capabilities, we 
would need an even more precise measurement of M. than 
possible with 300 low-z SNIa. With such a measurement, 
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the determination of a would be very accurate indeed. 

One can argue that the inverse power-law potential is 
not the most favoured using current data because it does 
not provide low enough values of wq. Sticking to inte- 
ger a one finds wq > —0.8 today, and values as high 
as a — 10 are already ruled out with present data. We 
have included them here to illustrate the method used 
and the possible accuracy of the SNAP satellite. In view 
of this, the most important aspect of this analysis is the 
ability to rule out the inverse power-law potential. As 
such we have shown that the method is effective and it 
is immediately generalizable to other potentials which 



mimic a cosmological constant for some choice of param- 
eter^). Our fitting procedure may also be used to rule 
out the inverse power-law potential in the case of other 
quintessence potentials. This is work in progress. 
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